Local Moment Formation in the Superconducting State of a Doped Mott Insulator 
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A microscopic theory is presented for the local moment formation near a non-magnetic impurity 
or a copper defect in high-T c superconductors. We use a renormalized meanfield theory of the t — J 
model for a doped Mott insulator and study the fully self-consistent, spatially unrestricted solutions 
of the d-wave superconducting (SC) state in both the spin S — and 5=1/2 sectors. We find 
a transition from the singlet d-wave SC state to a spin doublet SC state when the renormalized 
exchange coupling exceeds a doping dependent critical value. The induced S — 1/2 moment is 
staggered and localized around the impurity. It arises from the binding of an S = 1/2 nodal 
quasiparticle excitation to the impurity. The local density of states spectrum is calculated and 
connections to NMR and STM experiments are discussed. 
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The local electronic structure near non-magnetic im- 
purities in high-T c superconductors has attracted wide 
attention recently. A series of NMR experiments in 
YBa2Cu307_j; (YBCO) have shown that with the sub- 
stitution of nonmagnetic Zn/Li for the Cu atoms in the 
Cu0 2 plane, an S = 1/2 staggered magnetic moment is 
generated on the Cu ions in the vicinity of the impurity 
sites |l|,^.||,^|J|,^,|^,|| . Below the superconducting tran- 
sition temperature, the moment is partially screened in 
optimally doped YBCO, but remains unscreened down 
to the lowest temperatures in the underdoped regime 
Low-temperature STM experiments |J have di- 
rectly observed the local electronic structure around the 
Zn impurity atoms on the surface of superconducting 
Bi 2 Sr 2 CaCu 2 8+2; (BSCCO). The tunneling differential 
conductance above the Zn site exhibits a sharp in-gap 
resonance peak near zero bias (-1.5meV) for electrons 
tunneling out of the sample || . 

While the presence of a resonance impurity state is 
consistent with the strong local potential scattering near 
the impurity 10 ll[] , the latter does not account for the 
local magnetic moment induced by nonmagnetic impuri- 
ties observed by NMR. This has led to attempts to iden- 
tify the low energy conductance peak with the Kondo 
resonance resulting from the Kondo screening of the lo- 
cal moment by the superconducting (SC) electrons Jl2| . 
However, the formation of the local moment itself near 
nonmagnetic impurities in a d-wave superconductor has 
not been addressed. In this paper, we present a theory 
which simultaneously explains the local moment forma- 
tion and predicts a strong low energy resonance in the 
local density of states (LDOS) spectrum. 

We study the local moment formation in d-wave su- 
perconductors in proximity to Mott insulators. In the 
resonance- valence-bond (RVB) picture p3| , the electrons 
are paired into spin singlets, which are mobilized by the 
doped holes and condense into the SC state below T c . 
A nonmagnetic impurity atom such as Zn/Li or a Cu 



vacancy locally breaks a singlet, creating an unpaired 
spin. This extra S — 1/2 moment can be self-consistently 
trapped by the local impurity potential, forming a local 
moment confined to the neighboring sites of the impu- 
rity. Our strategy is as follows. The d-wave BCS ground 
state |BCS) is a spin singlet and the low-lying excitations 
are the nodal quasiparticles. A single quasiparticle exci- 
tation, 7t|BCS), carries a spin S — 1/2. The unpaired 
spin induced by the impurity corresponds to a superposi- 
tion of such one-7 excitations. We refer to such a doublet 
state as a 7BCS state. If the self-consistently determined 
energy of the 7BCS state is lower than that of the BCS 
ground state, the true ground state is a doublet with a 
net S — 1/2 moment. Below, we show that the physi- 
cal picture described above can be realized by a renor- 
malized slave boson meanfield theory of the t — J model. 
We carry out fully self-consistent, spatially unrestricted 
calculations of the valence bond, SC order parameter, 
local doping concentration (LDC), and the spin density 
distribution in the doublet 7BCS state and compare its 
energy to that of the singlet d-wave BCS state. We find 
that binding occurs and the doublet 7BCS state becomes 
the true ground state when the renormalized exchange 
coupling a J exceeds a nonzero doping dependent critical 
value. The S = 1/2 moment exhibits a staggered pattern 
localized around the impurity in agreement with the re- 
sults derived from the NMR experiments. We determine 
the phase boundary in the aJ versus doping plane, ob- 
tain the tunneling density of states in the doublet 7BCS 
state and compare to the results of STM experiments. 

We emphasize the importance for studying ground 
states with good spin quantum numbers. In a recent 
work pi , Gutzwiller approximation of the t-J model 
was used to calculate the induced moment near a non- 
magnetic impurity. However, the authors did not seem 
to have considered states with well defined spin quantum 
numbers. They reported a net moment of a fraction of an 
electron spin which was obtained due to a partial polar- 
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ization of the singlet BCS state in the meanfield theory. 
Therefore, the binding of an spin-1/2 excitation to the 
impurity i.e. the formation of the local moment has not 
been addressed. 

We begin with the t — J model on a square lattice in 
the presence of the long-range Coulomb interaction, 

H = ~ (tijCivCj" + h - c -) + J E ( S * ' S i~ 4 n * n j) 
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Here c ia creates an electron and the spin operator S; = 
{\)c ia a a pCij3. The sum over i,j in the hopping term in- 
cludes the nearest neighbors (n.n.) where Uj = t and 
the next n.n. where tij = t' . To model a nonmagnetic 
impurity or a Cu defect, we remove a site on the L x L 
square lattice as indicated by the primes on the sum- 
mations in Eq. (|l|). The total number of sites is then 
N s = L X L — 1. The long-range Coulomb interaction 
strength V c in Eq. (Q) is given by V c — e 2 /47rea ~ t 
p5[ . Eq. (Q) is only physical in the limit of strong on- 
site Coulomb repulsion, i.e. under the constraint of no 
double occupancy, »j = ct^c^ < 1. Such a projected 
Hilbert space can be treated in the slave-boson formula- 
tion by writing c[ a — fj a bi, where is a spin-carrying 
fermion and hi a spinless boson fill with the constraints 
fiafia + o[bi = 1. The AF spin-exchange term is decou- 
pled according to fL7|] 
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where Ay and Xij are the pairing and valence bond order 
parameters on each n.n. bond respectively, 
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In the meanfield SC state, the bosons locally condense 
[Ea,16| into the coherent state, 
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The LDC is given by x, = {^ b \b\bi\^ b ) = 
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The 



constraint 



= 1 — m is enforced on average at every 
site by locally adjusting the chemical potential. 

It is known that fluctuations beyond such a meanfield 
RVB state will mediate a residual direct spin-spin inter- 
action, which can be approximately accounted for phe- 
nomenalogically by adding to the meanfield theory an 
interaction term, 



H a = aJ^2 ( S i) ' S j; 
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where a is a phenomenalogical parameter between zero 
and one. This renormalized meanfield approach was used 
recently for a successful description of the spin dynamics 
in YBCO |l7fl . The fully self-consistent, spatially unre- 
stricted solution will include nonlinear magnetic effects 
induced by H a near the impurity. 

We diagonalize the theory in real space by writing 
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where 7^, n = 1,...,2N S creates a Bogoliubov quasi- 
particle with energy E n and wave- function (u n i,v n i) as 
solutions of the Bogoliubov-de Gennes (BdG) equations. 
Since the particle-hole symmetry in the quasiparticle ex- 
citation spectrum, which reflects the degeneracy in spin 
excitations, is broken in the presence of a net spin mo- 
ment, we keep all independent eigenstates with E n < 
and E n > and construct the ground state according to 
the occupation of the fermionic quasiparticle orbitals: 

Tt occ Tt occ N s 

\* G ,n occ ) = n 7ilo) = n D<4 + <ifnM, (6) 



where n occ labels the last 7-quasiparticle level to be 
occupied. The spin-singlet d-wave BCS state corre- 
sponds to |BCS) = I^Gj^occ = N s ), whereas the 7BCS 
state is given by | 7 BCS> = \^ G ,n occ = N s + 1), 
which is a global doublet with a net S = 1/2 mo- 
ment. These states must be determined by minimiz- 
ing (*G,"occ|^mf + -ff Q |*G,«occ), leading to the self- 
consistency equations 
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The LDC and the local spin density are given by 
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For l^c^occ) to be a self-consistent ground state, it re- 
quires a lowest total energy and, in addition, E Uacc < 0. 

The BdG equation and the spatially unrestricted self- 
consistency equations can be solved by numerical itera- 
tions starting from a random set of initial conditions. We 
use J as the unit of energy and set t/J — 3, V c /J = 5. 
The impurity/defect is located at (L/2, L/2). We first 
consider t' = and focus on the nature of the ground 
state. The energies per site for the BCS and the 7BCS 
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states are shown in Fig. 1 as a function of the renormal- 
ized spin-spin interaction aJ for L — 24 at an average 
doping x — 0.14. While that of the singlet state is in- 
dependent of aJ, the energy of the doublet 7BCS-state 
is reduced with increasing aJ and drops abruptly below 
that of the singlet state near a critical value a c , indicating 
that the doublet state can gain substantial spin-exchange 
energy through the formation of the magnetic moment. 
Notice that although the doublet state energy is already 
lower than that of the singlet state for a less than a c , it 
is not a self-consistent ground state unless the last occu- 
pied N s + 1 quasiparticle level is pulled below the Fermi 
level (see Eq. (||)). To determine the value of a c at which 
a phase transition into the 7BCS-state occurs, we calcu- 
late the energy of the n occ = N s + 1-th 7-quasiparticle 
level E noca as a function of a J and identify a c ~ 0.68 at 
x = 0.14 by the point where it crosses the Fermi energy. 
The phase diagram in the plane of a J versus doping x is 
plotted in the inset of Fig. 1, which shows that the forma- 
tion of the local moment requires a smaller renormalized 
spin interaction in the underdoped regime. 
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FIG. 1. Energy per site of the singlet and the doublet 
states as a function of a J on a 24 x 24 lattice at x = 0.14. 
The critical value of the exchange coupling is a c ~ 0.68. 
Inset: phase diagram in the aJ versus doping x plane. 

Next, we turn to the spatial distribution of the im- 
purity/defect induced local moment in the 7BCS-state. 
We choose large systems with L — 32 and set t' = —t/3 
and a = 0.75 > a c at x — 0.14. The self-consistent 
solutions of the LDC Xi and the d-wave gap amplitude, 
Af = (Ai,i+i + A M _ £ - A i)i+S - A M _ 5 )/4, are shown in 
Figs. 2a and 2b as 2D maps. Both Xi and Af are strongly 
suppressed in the vicinity of the impurity/defect, but 
quickly resume their bulk values away from it. The va- 
lence bond Xij shows only moderate oscillations near the 
impurity with no signatures for the emergence of a stag- 
gered flux or current. The spatial distribution of the mag- 
netic moment Sf , having a net magnitude J^. Sf = 1/2, 
is plotted in Fig. 2c. It shows a staggered pattern local- 



ized around the impurity in good agreement with that 
derived from the NMR experiments ffil. The majority 
of the net moment resides on the four nearest neighbors 
of the impurity. Interestingly, along the nodal directions 
(diagonals) of the d-wave gap, weak incommensurate spin 
density wave (SDW) oscillations arise with a periodicity 
of approximately 8 unit cells. This feature is more clearly 
seen in the 3D plot of the staggered magnetic moment 
shown in Fig. 2d. The dependence of the SDW period- 
icity on the value of t' and the doping x is consistent 
with the behavior of the bulk spin susceptibility which 
peaks around similar incommensurate wave- vectors [[l7| . 
Physically, an impurity /defect induces a local moment 
in its immediate vicinity, which subsequently leaks into 
the bulk via the short-range spin-spin correlation in the 
short-range RVB state. 




FIG.°2. Serf-conS9stcn1 3< solution of the 7BCS-?ftate 1 with 
a net S = 1/2 local moment on a 32 x 32 lattice with 
a = 0.75 > a c and t' = —t/3. The impurity /defect is lo- 
cated at (15, 15). (a) LDC xt. (b) d-wave gap amplitude 
Af. (c) Distribution of the local moment S*. (d) 3D-plot 
of the staggered magnetic moment. 

We next study the LDOS spectrum measured by STM. 
Note that the net spin polarization in the doublet ground 
state flips between S z = ±1/2, which should be identified 
with the two degenerate states of an S = 1/2 quantum 
impurity. We set S z = 1/2 and do not consider here the 
effect of Kondo coupling between such an extended local 
moment and the quasiparticles which can be quite dif- 
ferent from the conventional Kondo problem in d-wave 
superconductors. The spectral weight of the electron is a 
product of those associated with the boson and fermion 
states given in Eqs. @ and @ respectively. We derive 
the LDOS for the up and down spin electrons, 
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Nufu) = [£ A+ \v ni \ 2 + A nl\ v ni\ 2 ]S(u + E n ), 

n n— n occ + l 

where A^ a is the boson coherent state overlap for adding 
(+) or removing (— ) an electron with spin-cr accom- 
panied by adding/removing a 7-quasiparticle from the 
n-th fermionic orbital. In order to maintain the oc- 
cupation constraint in the final states, the correspond- 
ing boson coherent states, defined as (j>f nlJ , must satisfy 

GO 2 = ^nl? = # ±J«ni| a T \Vni\ 2 . Using Eq. §, 
we obtain the spectral weight for the bosons, 
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FIG. 3. The LDOS spectra at the n.n. sites of the 
impurity /defect (solid lines) and near the corner (dashed 
lines) of 32 x 32 lattices at x = 0.14. (a) t' = -0.45t and 
aj = 0.7, close to the phase boundary, (b) t' = — 1/3 and 
aj = 0.75, the same system as in Fig. 2. Insets: details 
of the resonance line-shape near zero bias with reduced 
thermal broadening. 

Due to the localized nature of the impurity states, the 
order parameters in Eqs. (|7|-|Tc|) are insensitive to the 
boundary conditions for large system sizes. We have cal- 
culated Ni(w) on 32 x 32 lattices and averaged the results 
over a set of different boundary conditions to reduce the 
finite-size effects. The complex boundary phases along 
the x and y directions are typically 2ir(i x /5,i y /5) where 
i x , i y = 0, ... ,4. The results are not sensitive to the de- 
tailed choice of the phases. In Fig. 3, we plot the LDOS 
as a function of the bias voltage for two ground states, 
one is close to the phase boundary (Fig. 3a) and the other 
(Fig. 3b) is the same system used in Fig. 2 in the local 
moment phase. On the n.n. site of the impurity/defect, 
the LDOS (solid lines) exhibits a near zero-bias resonance 
with its peak located on the occupied side of the spec- 
trum. Comparing to the LDOS obtained near the corner 



of the system (dashed lines) , the coherence peaks at the 
gap edges have been strongly suppressed around the im- 
purity/defect. These results are in good agreement with 
STM observations The LDOS spectra with reduced 
thermal broadening are plotted in the insets of Fig. 3 
to reveal the details of the resonance line-shape. The 
inset in Fig. 3a is in good agreement with current experi- 
ments |9). However, comparison with the inset in Fig. 3b 
shows that the resonance line-shape is sensitive to the 
parameters and show several peaks asymmetrically dis- 
tributed around zero bias. The low energy spectra can 
be qualitatively understood by considering the last two 
occupied quasiparticle levels: £jv s and £^jv s +i- The reso- 
nance state at these two energies is occupied in the local 
moment phase which should in principle lead to four reso- 
nance peaks centered at ±En s and ±En s+ i in the LDOS 
spectrum through their particle-hole reflections. How- 
ever, since the SC coherence is significantly reduced near 
the impurity /defect, the spectral weights have the prop- 
erty |ujvj 2 <C \vn s \ 2 whereas |-uat s+ i| 2 3> |vat s+ i| 2 lead- 
ing to two predominant peaks near zero bias. The peak at 
the negative energy £W s +i (positive energy \En b |) corre- 
sponds to removing (adding) an up (down) spin electron 
from the resonance state. 

In conclusion, we have presented a microscopic theory 
that describes on equal footing the formation of local mo- 
ment observed by NMR and the local electronic structure 
measured by STM near nonmagnetic impurities in high- 
T c superconductors. The existence of the low energy 
resonance is a robust feature of our model, even though 
the detailed structure of the resonance may vary. Unless 
the surface of BSCCO is very close to the onset of local 
moment formation such that the detailed resonance peak 
structure becomes difficult to observe, our model predicts 
multiple peak structures for better resolved resonances in 
the local moment phase. Moreover, since recent STM re- 
vealed a spatially inhomogeneous electronic structure in 
BSCCO [[l8| that can be interpreted as off-stoichiometry 
carrier doping induced spatial variations in the LDC [Gl| , 
the present theory would predict a distribution of the res- 
onance structure over different Zn-impurity sites. 
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